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摘 要 : 为 了 解决 粒子 滤波 多 说 话 人 跟踪 过 程 中 粒子 易 发 散 导 致 多 目标 跟踪 精度 低 的 问题 ， 提 出 了 并 行 粒子 滤波 和 基 
于 GPU 的 K- 均 值 聚 类 的 多 声 源 定位 方法 。 该 方法 首先 分 析 了 粒子 滤波 在 实现 多 目标 跟踪 时 ， 进 行 数据 关联 的 过 程 产 
生 较 大 的 计算 量 ， 并 且 出 现 多 个 目标 时 ， 粒 子 会 逐渐 发 散 。 针 对 计算 量 大 和 粒子 发 散 的 问题 ， 提 出 了 一 种 并 行 粒子 滤 
波 和 开 - 均 值 聚 类 的 方法 。 实 验 表 明 ， 随 着 粒子 数 和 目标 数 的 增加 ， 计 算 量 以 指数 增加 ， 并 且 粒 子 发 散 严重 ， 采 用 基于 
GPU 的 开 - 均 值 聚 类 方法 的 粒子 滤波 多 说 话 人 跟踪 方法 ， 相 比 传 统 粒子 滤波 跟踪 方法 具有 更 收敛 的 粒子 集 并 且 跟 踪 精 
度 较 高 。 
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Particle filter multi-speakers tracking algorithm based on GPU and its application 
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Abstract: In order to solve the problem of low accuracy of multi-target tracking, particles in particle filter are easy to disperse 


in the process of multi-speaker tracking. This paper presented a parallel particle filter algorithm and GPU-based K-means 


clustering multi-source localization method. The method first analyzed the particle filter to achieve multi-target tracking, data 


association process have a large amount of computation, and the particles gradually diverge with the emergence of multiple 


targets. In order to solve the problem of large amount of computation and particle divergence, this paper proposed a method of 


parallel particle filter and k-means clustering. Experiments show that, as the number of particles and the number of targets 


increases, the amount of computation increases exponentially and the particles scatter seriously. In this paper, by using the GPU- 


based k-means clustering method, the particle filter multi-speaker tracking method has more convergent particle sets and higher 


tracking accuracy than the traditional particle filter tracking method. 
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粒子 滤波 方法 09 非 高 斯 非 线性 的 系统 模型 ， 这 更 加 接近 真是 场 


0 引言 景 ， 并 将 时 延 估 计 作为 观测 量 ， 滤 波 后 直接 得 到 声 源 的 位 置 ， 

说 话 人 运动 过 程 0 习 具有 较 强 的 随机 性 , 通过 采用 合适 的 模 该 方法 解决 了 非 线性 非 高 斯 模型 下 的 系统 状态 估计 问题 。 实 际 
型 来 描述 说 话 人 的 运动 过 程 ， 有 研究 人 员 提 出 用 动态 空间 模型 场景 通常 是 非 线 性 非 高 斯 的 ， 因 此 粒子 滤波 在 说 话 人 的 定位 与 
的 方法 描述 说 话 人 的 运动 , 并 利用 滤波 理论 B4 解 决 说 话 人 跟踪 跟踪 应 用 最 广泛 。 后 续 研 究 中 ， 许 多 研究 者 针对 粒子 滤波 在 说 
的 问题 ， 动 态 空间 模型 提供 了 说 话 人 运动 规律 的 先 验 信 息 ， 利 话 人 跟踪 的 缺陷 提出 改进 ， 改 进 方向 有 三 个 ， 一 是 定位 特征 的 
用 先 验 信息 有 利于 提高 说 话 人 定位 的 精度 。 在 说 话 人 的 定位 与 提取 [由 二 是 粒子 滤波 算法 的 改进 433， 针对 算法 的 不 足 提出 改 
跟踪 中 , 卡尔 曼 滤 波 方法 569 解决 了 高 斯 线性 模型 下 的 系统 状态 进 并 应 用 到 说 话 人 跟踪 中 ; 三 是 采用 更 合适 的 动态 空间 模型 来 
估计 问题 ; 扩展 卡尔 曼 滤 波 方法 [9 时 延 估 计 作 为 观测 量 进行 非 描述 多 说 话 人 运动 过 程 。 另 外 传统 粒子 滤波 在 多 目标 跟踪 时 ， 
线性 滤波 ， 解 决 了 高 斯 非 线 性 模型 下 的 系统 状态 估计 问题 ， 无 通过 使 用 多 个 粒子 滤波 器 ， 每 个 滤波 器 跟踪 一 个 目标 ， 利 用 数 
际 卡尔 曼 滤波 方法 四 解决 了 高 斯 非 线性 模型 的 状态 估计 问题 。 据 关联 方法 将 不 同 的 量 测 值 与 目标 进行 相对 应 ， 从 而 获得 每 个 
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目标 的 特征 ， 每 个 目标 的 跟踪 是 一 个 独立 的 过 程 。 关 联 过 程 会 
产生 较 大 的 计算 量 ， 并 且 粒 子 滤 波 在 多 目标 跟踪 中 易 产生 粒子 
发 散 和 跟踪 精度 低 的 问题 。 

针对 计算 量 大 的 问题 ， 本 文 引 入 了 并 行 的 思想 ， 将 粒子 滤 
波 并 行 实现 ， 并 且 针 对 粒子 滤波 进行 多 目标 跟踪 过 程 中 ， 重 采 
样 后 的 粒子 易 发 散 的 问题 ， 本 文采 用 基于 GPU 的 -均值 聚 类 
方法 3-161|， 将 K- 均 值 率 类 算法 中 计算 密集 的 部 分 交 给 GPU 端 
执行 计算 。 
1 ”说 话 人 运动 模型 

多 说 话 人 运动 模型 用 于 表示 移动 的 说 话 人 在 智能 会 议 室 环 
境 里 位 置信 息 的 变化 ， 其 中 广泛 采用 的 模型 为 Langevin 模型 。 
该 模型 综合 了 各 类 随机 运动 的 特点 ， 在 声 源 的 定位 与 跟踪 问题 
中 ， 该 模型 能 够 较 好 地 描述 说 话 人 的 运动 过 程 ， 因 此 被 广泛 采 


pe 


声 源 在 半 轴 方向 的 运动 方程 如 下 


,=[X, yd,] (1) 
T=ax thr 
X=X, +ATX, 
a O) 
a.=e 


b=v, Vl-a? 
其 中 : x 和 这 分 别 表 示 马 轴 方 向 的 位 置 和 速度 ，y 和 分别 表 
示 y 轴 方向 的 位 F， 是 正 态 分 布 的 随机 变量 


和 速度 。 


AT=D 
大 


常数 ， 为 稳 态 均 方 根 速度 ， 本 文中 8.=104s ，v =101ms 。 


为 帧 间 的 时 间 间 隔 , 工 是 帧 长 ， 太 是 采用 频率 ，A 为 


2 多 目标 模型 


根据 目标 运动 状态 和 数量 的 变化 , 目标 的 运动 模型 有 四 种 : 
目标 运动 ， 目 标 消 失 ， 目 标 产 生 和 目标 衍生 ， 如 图 1 所 示 。 


时 刻 上 时 刻 k+1 
Pe 


状态 从 xza 到 xs 的 概 府 


状态 为 xia 的 目标 消失 的 
概率 


状态 为 xi 的 目标 产生 的 
随机 集 及 了 


产生 的 随机 集 /7 


标的 运动 模型 


2.1 基于 粒子 滤波 和 K- 均 值 聚 类 的 多 说 话 人 跟踪 方法 


传统 粒子 滤波 在 进行 多 目标 跟踪 时 ， 通 过 使 用 多 个 粒子 滤 
波 器 ， 每 个 滤波 器 跟踪 一 个 目标 ， 利 用 数据 关联 方法 将 不 同 的 


量 测 值 与 目标 进行 相对 应 ， 从 而 获得 每 个 目标 的 特征 ， 每 个 
标的 跟踪 是 一 个 独立 的 过 程 。 关 联 过 程 会 产生 较 大 的 计算 量 ， 
上 且 传 统 粒子 滤波 在 对 多 个 目标 进行 跟踪 时 ， 随 着 目标 数量 的 
增多 ， 粒 子 易 出 现 发 散 的 现象 。 针 对 粒子 发 散 和 计算 量 大 的 问 
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题 ， 本 文 引 入 并 行 策略 ， 将 粒子 滤波 与 KK- 均值 聚 类 算法 并 行 实 
现 ， 从 而 减少 计算 复杂 度 ， 从 而 达到 提高 粒子 滤波 跟踪 算法 跟 
踪 精 度 的 目标 。 
2.2 基于 GPU 的 K- 均 值 聚 类 算法 

K -均值 聚 类 算法 核心 思想 是 : 首先 从 所 给 的 n 个 对 象 中 选 
取 大 个 对 象 作为 初始 聚 类 中 心 点 ， 对 剩余 的 其 他 对 象 计算 它们 
与 所 选 的 大 个 中 心 点 的 相似 度 ， 然 后 再 重新 计算 所 获 聚 类 的 聚 
类 中 心 。 算 法 基本 流程 如 下 : 

a) 从 粒子 集 X 中 随机 选择 大 个 对 象 作为 初始 机 
Cl,C2s°*'s CL ° 

b) 将 对 象 xG=1,2,3,…,n) 按照 距离 分 配给 最 近 的 一 个 聚 
类 中 心 。x 和 聚 类 中 心 的 距离 定义 为 


lx -ol | -ey G) 


其 中 : m 为 数据 的 个 数 。 
c) 重新 计算 每 个 簇 中 新 的 聚 类 中 心 c， 
4 , j=1,2,.…,K (4) 


j xes; 


其 中 N) 为 第 j 个 聚 类 艇 9 中 对 象 的 个 数 。 
d) 重复 b) c) 直到 大 个 聚 类 中 心 不 再 变化 ， 即 准则 函 


二 


聚 类 中 心 


数 收 


基于 GPU 的 天- 均值 主要 思想 是 数据 并 行 ， 将 传统 天- 均值 
聚 类 算法 中 计算 密集 的 部 分 放 在 设备 端 进行 计算 ， 从 而 实现 提 
高 性 能 的 目的 。 数据 对 象 分 配 和 -均值 中 心 点 的 计算 执行 多 次 ， 
且 不 同 的 数据 独立 ,将 数据 对 象 分 配 和 -均值 中 心 点 计算 可 
以 分 离 成 两 个 内 核 函数 ， 该 函数 由 大 量 线程 组 成 ， 在 主机 端 调 
用 ， 在 设备 端 执行 。 在 任务 分 配 时 ， 主 机 负责 将 k 对 象 放置 到 
被 聚 类 的 对 象 所 表示 的 空间 ， 并 重新 安排 所 有 的 数据 对 象 与 控 
制 和 迭代 过 程 ， 而 设备 端 进行 数据 并 行 密 集 计 算 。 在 数据 对 象 赋 
值 之 后 ， 得 到 每 个 数据 点 的 艇 标签 ， 用 于 计算 一 个 聚 类 中 心 。 
本 文 将 聚 类 标签 从 设备 端 下 载 到 主机 端 ， 主 机 端 重 排 所 有 的 数 
据 并 计算 每 个 群集 所 包含 的 数据 对 象 的 数目 。 

算法 1 基于 GPU 的 -均值 聚 类 算法 

1 if threadId=0 then 


Im 
Dn 


20, c+ 

3 endif 

4 Synchronize threads 

5 repeat 

6 forall GC {C+ do 
7 C, Cr + 

8 end for 

9 synchronize threads 


10if threadld=0 then 


11 forall GC, C+ do 
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14 forall Cj cdo 


1 
15 C; <€——=—C, 


16 end for 

17 if convergence then 

18 signal threads to terminate 
19 end if 


20 until convergence 


在 任务 分 配 时 , 主机 负责 将 K 对 象 放置 到 被 聚 类 的 对 象 所 


表示 的 空间 中 ， 并 重新 安排 所 有 的 数据 对 象 与 控制 迭代 过 程 ， 
而 设备 端 进行 数据 并 行 密集 计算 。 在 数据 对 象 赋值 之 后 ， 得 到 


每 个 数据 点 的 艇 标签， 用 于 计算 一 个 聚 类 中 心 。 最 简单 的 思路 
是 读 取 所 有 的 数据 对 象 和 确定 数据 点 的 聚 类 中 心 。 然 而 ， 大 量 
的 条 件 语句 不 适用 于 GPU 处 理 器 模型 。 因 此 本 文 增加 一 个 过 
程 ， 这 个 过 程 包括 将 聚 类 标签 从 设备 端 下 载 到 主机 端 ， 主 机 端 
重 排 所 有 的 数据 并 计算 每 个 群集 所 包含 的 数据 对 象 的 数目 。 在 
数据 存储 时 ， 数 据 对 象 和 中 心 点 以 动态 数组 的 方式 在 设备 端 储 
存 。 本 文 将 K- 均 值 聚 类 算法 的 所 有 参数 存储 在 在 全 局 内 存 中 ， 
因为 常数 存储 器 和 纹理 存储 器 均 为 只 读 的 存储 器 ， 且 存储 内 存 
只 有 64kb， 内 存 太 小 以 致 于 不 够 容纳 K 均值 聚 类 的 所 有 参数 。 
于 GPU 是 一 种 共享 内 存 的 处 理 器 架构 ， 处 理 器 被 称 为 
线程 ， RE 每 个 线程 都 被 限定 在 0 至 tl 之 间 

是 线程 总 数 。 通 常情 况 下 ， 线 程 0 被 认为 是 主线 程 ， 所 有 划 
他 的 线程 被 认为 是 从 线程 (任何 一 个 线程 都 可 以 被 认为 是 主线 
程 ， 其 他 为 从 线程 )。 线 程 共享 数据 点 集 所 在 的 一 部 分 内 存 ， 同 
样 还 包括 当前 的 中 心 点 集 w ， 还 有 计算 所 有 内 部 中 心 点 直接 的 
距离 n 矩阵 A 和 对 该 矩阵 A 进行 排序 ， 获 得 下 标 排列 矩阵 B。 
每 一 个 线程 还 附加 了 三 种 内 存 ， 供 内 部 数据 使 用 。 假 设 对 封闭 
的 架构 访问 是 可 行 的 。 映 射 到 GPU 的 编程 描述 如 下 : 首先 主线 
程 和 串 行 算法 一 样 初始 化 所 有 的 中 心 点 ， 然 后 n 被 分 割 成 子 集 
n。 这 仅仅 是 标记 阶段 ， 每 个 线程 执行 赋予 n 单独 的 处 理 器 的 
坐标 和 计算 范围 。 所 有 线程 执行 标记 阶段 都 是 针对 于 它们 各 自 
的 n 划分 。 每 个 数据 点 的 标记 都 被 存储 在 一 个 n 维 向 量 的 组 
件 7 中 。 这 样 便 消除 了 当 更 新 聚 类 和 简化 记录 过 程 中 ， 同 时 写 
的 问题 
在 标记 阶段 之 后 ， 线 程 间 需要 进行 同步 操作 ， 以 确保 所 有 
中 心 点 更 新 的 阶段 数据 已 经 是 可 以 使 用 的 。 之 后 ， 中 心 点 更 新 
阶段 可 以 通过 规约 操作 来 执行 。 

然而 ,为 了 简化 ,假设 主线 程 执行 这 个 阶段 是 顺序 执行 的 。 
主线 程 不 是 欠 代 全 部 中 心 点 ， 而 是 迭代 部 分 计算 的 新 中 心 点 的 
全 部 标记 。 mj 维 向 量 mj 在 一 次 迭代 中 被 更 新 ， 向 量 的 每 
pe ee 案 类 C 的 数值 。 执 行 下 一 


个 循环 遍历 所 有 中 心 点 ， 用 j, 处 理 C, 来 获得 得 最 终 的 中 心 点 。 


收敛 判断 同样 是 由 主线 程 决 定 ， 主 线程 通过 检查 最 好 一 个 标记 


ChinaXiv 合 作 期 十 


等 : en da 算法 及 其 应 用 


划 ”法 
日 二， 


阶段 在 聚 类 算法 中 是 否 有 改变 来 确定 。 只 要 达到 收敛 的 条 件 ， 
从 线程 由 主线 程 来 停滞 其 执行 。 
基于 GPU 的 并 行 K- 均 值 聚 类 算法 流程 如 图 2 所 示 。 


聚 类 标签 下 载 


k 均 什 中 心 点 下 载 


图 2 基于 GPU 的 KK- 均值 算法 


2.3 定位 函数 
本 文 在 提取 定位 特征 时 ， 利 用 环形 阵列 ， 并 结合 可 控 功 率 
响应 方法 提取 多 说 话 人 的 方位 角 ， 即 


P(0)=Y > | yOX OX er do G5) 


通过 放置 两 个 环形 麦克 风 阵 列 ， 并 利用 这 两 个 环形 阵列 来 
提取 多 说 话 人 的 定位 特征 “图 3)。 


3 双环 形 阵列 示意 图 


设 坐 标的 原点 O 位 于 两 个 环形 阵列 的 中 间 位 置 ， 两 个 环形 


阵列 的 中 心 坐标 分 别 为 6 -Co wy ) 和 局 =(xo ye )， 声 源 的 位 置 


华 标 为 元 =(x,y) ， 则 声 源 对 应 于 两 个 环形 阵列 的 方位 角 8,6, 分 
别 满足 : 
XXo 
Crxo) +O-y0) 2 
X-Xo, 
a We 
2.4 基于 粒子 滤波 多 说 话 人 跟踪 算法 
本 节 利 用 并 行 粒子 滤波 算法 对 多 说 话 人 的 状态 进行 跟踪 ， 
其 中 的 运动 模型 采用 Langevin 模型 , 似 然 函数 采用 高 斯 似 然 函 
数 , 将 并 行 粒子 滤波 与 并 行 K- 均 值 聚 类 方法 相 结 合 实现 多 声 源 
定位 跟踪 的 步骤 如 下 : 
a) 利用 初始 分 布 采样 粒子 {@i,i=1,…， 


Q=arccos( 


=arccos( 


N} ， 对 应 权 值 为 
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{ws =1/N,i=1,.…, N}。 
b) 根据 Langevin 运动 模型 预测 最 新 时 刻 粒子 的 位 置 分 布 。 
c) 估计 每 个 粒子 的 方位 角 计算 似 然 函数 。 方 位 角 根 据 式 
(6)0) 计 算 ， 然 后 计算 对 应 的 似 然 函数 。 


d) 将 似 然 函数 计算 的 结果 作为 粒子 的 权 值 ， 并 对 权 值 进行 
归 一 化 处 理 ， 
~ we 
Wi wp (8) 
利用 归 一 化 后 的 权 值 对 粒子 进行 重 采样 。 


e) 针对 多 说 话 人 跟踪 问题 ， 需 要 利用 并 行 玉 - 均 值 聚 类 来 
确定 哪些 粒子 来 自 于 同一 个 跟踪 目标 ， 在 并 行 玉 -均值 聚 类 时 ， 
说 话 人 的 数目 己 知 。 在 聚 类 后 说 话 人 的 位 置 估计 为 


划 法 
日 二 ， 


其 中 : &, 是 第 的 位 
NN, 为 粒子 数 


3 ”实验 结果 与 分 析 


实验 平台 中 硬件 ; 


Visual Studio 2013， 实 验 


等 : 基于 GPU ee 
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QW 为 对 应 于 第 p 个 说 话 人 的 粒子 ， 


台式 机 和 GTX960 显卡 ， 软 件 平 台 为 
伟人 硬件 相关 参数 如 表 2 所 示 。 


表 2 实验 硬件 相关 参数 
GPU CPU 
GPU GTX960 CPU i5-4460 
核 数 1024 核 数 4/8 
时 钟 率 ”1.3 GHz ”时 钟 率 ”3.2 GHz 


在 粒子 滤波 的 多 声 源 跟踪 


中 引入 


是 对 粒子 滤波 重 采样 后 的 粒子 作 午 
的 粒子 归 为 一 类 , 因此 实验 首先 给 出 了 


行 K- 均 值 聚 类 的 目的 
聚 类 分 析 ， 将 属于 同一 个 目标 
行 玉 -均值 聚 类 算法 对 


CQ = -Do (9) 
粒子 的 聚 类 过 程 ， 如 图 5 所 示 。 
5 5 
六 | | xu 
: ~ 4 “i 
3 3l 
村 2 革 
s 2 E34 
六 六 于 理 
0 0 
5 0 1 3 4 5 0 1 2 3 4 5 
x 轴 X 轴 
(al) 重 采样 前 的 粒子 分 布 (bl) 重 采样 后 粒子 分 布 (cl1)K- 均 值 聚 类 后 的 结 
5 5 5 
4. 4 4 i 
3 3r 3 上 
四 想 四 
> 中 2. FE 于 
1 1 -i 1 本 
00 1 2 3 4 5 00 1 2 3 4 5 "0 1 2 3 4 5 
x 轴 x 四 
(a2) 重 采样 前 的 粒子 分 布 (b2) 重 采样 后 粒子 分 布 (c2) 并 行 K- 均 值 聚 类 后 的 结 
图 5 并 行 执行 时 不 同 阶段 粒子 分 布 的 变化 


从 图 6 中 可 知 ， 在 重 采样 前 粒子 分 布 相对 分 散 ， 而 在 重 采 
样 后 一 些 权 值 低 的 粒子 被 大 权 值 的 粒子 蔡 换 ， 重 采样 后 的 粒子 
分 布 相对 比较 集中 ， 虽 然 从 粒子 的 分 布 图 中 可 以 明显 地 观察 到 
每 个 粒子 所 属 的 目标 ， 但 是 粒子 滤波 算法 本 身 无 法 对 这 些 粒子 
归 类 ,图 6(c) 则 给 出 了 并 行 K- 均 值 聚 类 算法 对 粒子 做 分 类 的 结 
果 , 显然 属于 不 同 目标 的 粒子 被 很 明显 地 区 分 开 , 相 比 K- 均 值 
聚 类 算法 ， 并 行 K- 均 值 聚 类 具有 更 明显 的 区 分 性 。 

为 验证 本 文 方法 在 智能 会 议 室 室内 场景 中 有 多 个 说 话 人 的 定位 
跟踪 效果 ， 实 验 场景 如 图 7 所 示 ， 智 能 会 议 室 的 噪声 及 混 响 设 
置 SVR = 204B,T =150ms ， 放 置 两 个 环形 麦克 风 阵 列 ， 房 间 中 
有 三 个 说 话 人 同时 发 声 ， 运 动 轨迹 分 别 为 直角 形 、 弧 形 和 直线 


型 ， 每 个 说 话 人 的 发 声 时 间 为 20 s。 


x 轴 


图 6 仿真 实验 的 环境 设置 图 
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(b) 并 行 粒子 滤波 + 聚 类 算法 的 跟踪 结果 


图 7 粒子 滤波 + 聚 类 算法 的 跟踪 结 


从 上 面 实 验 可 
人 跟踪 在 声 源 数 


得 出 , 并 行 粒 子 滤波 和 
的 有 效 数 目的 性 能 比 粒 子 滤 ; 


行 K- 均 值 的 多 说 话 


类 的 多 说 话 人 跟踪 精度 高 。 并 且 本 文 算法 在 高 强度 噪声 环境 下 ， 


不 仅 能 实时 确定 声 源 


本 效 数 


,也 能 准确 跟踪 多 说 话 人 的 位 置 。 


从 图 7 中 可 
下 ， 本 文 方法 可 以 跟踪 到 多 个 说 话 人 的 位 置 ， 在 说 话 人 数 


以 看 出 ， 在 多 个 说 话 人 轨迹 交叉 不 严重 的 情况 


率 
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环 


0 5 10 15 20 
时 间 (s) 
图 9 YY 轴 的 跟踪 结果 


表 3 两 种 算法 在 X 方 向 的 跟踪 误差 比较 
算法 蓝 绿 红 


PF+K 均值 跟踪 误差 率 (%) 0.9015 1.5046 13.112 
本 文 算法 跟踪 误差 率 (%) 0.8052 0.7493 5.475 


从 图 8、9 中 分 别 选 出 前 10 个 时 刻 的 X 轴 的 位 置 ， 作 误差 
比较 ， 即 用 跟踪 得 到 的 位 置 与 说 话 人 的 真实 位 值 取 相对 误差 


值 
表 
统 
跟 


的 


K- 均 值 聚 类 的 多 说 话 人 跟踪 效果 好 。 


， 分 别 用 图 8 中 的 蓝 绿 红 三 条 轨迹 作 误差 比较 ， 得 出 结果 如 
3 所 示 。 根 据 表 3 中 误差 率 可 知 ， 本 文 算法 的 误差 率 都 比 传 
粒子 滤波 和 天- 均值 聚 类 的 跟踪 算法 误差 率 低 , 即 本 文 算法 的 
踪 精 度 高 于 传统 跟踪 算法 ， 本 文 算 法 性 能 更 好 。 

综 上 所 述 ， 并 行 粒子 滤波 和 基于 GPU 的 K- 均 值 聚 类 算法 
多 说 话 人 跟踪 在 声 源 数 目的 有 效 数 目的 性 能 比 粒子 滤波 和 
且 本 文 算法 在 高 强度 品 


8 和 9 分 别 给 出 了 x 轴 和 y 轴 的 跟踪 轨迹 。 


x 轴 


0 5 10 15 20 
时 间 (s) 


图 8 X 轴 的 跟踪 结 


知 的 情况 下 , 并 行 K- 均 值 从 类 算法 可 以 更 有 效 地 对 重 采 样 粒 于 
进行 划分 ， 从 而 得 到 了 精度 较 高 的 多 说 话 人 位 置 估计 结果 。 图 


性 


声 环境 下 ， 不 仅 能 实时 确定 声 源 有 效 数 目 ， 也 能 准确 跟踪 多 说 
话 人 的 位 置 。 
和 区- 均值 聚 
的 值守 。 4 。 结束语 
通过 分 析 粒 子 滤波 在 多 说 话 人 跟踪 中 存在 的 问题 ， 本 文通 
过 引入 并 行 策略 ,将 粒子 滤波 算法 并 行 实现 , 降低 算法 复杂 度 。 
目 己 分 析 了 多 个 目标 出 现在 跟踪 区 域 时 , 粒子 出 现 的 发 散 的 问题 ， 
采用 基于 GPU 的 KK- 均值 聚 类 的 方法 来 解决 。 实 验 结 果 表 明 ， 
本 文 算法 有 效 提高 算法 在 GPU 中 硬件 资源 利用 率 ， 粒 子 收敛 


更 强 ， 本 文 算法 最 终 的 多 目标 跟踪 中 跟踪 精度 也 比 原始 粒子 


滤波 跟踪 算法 高 。 
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